*Replication Odesa Paper 20171125*

log using "Replication Odesa Paper 20171125", replace

use "/Users/olaonuch/Desktop/UCEPS_ODESA_20171125.dta", clear

*Variable Creation for Odesa Paper


*DV: WHO DID THE KILLING IN ODESA?*
*1 = pro-russians-did-it, 2 = pro-Ukrainians did it, 3 = responses that are neutral with respect to the first two narratives. V54othcod is a hand-coding by Henry Hale and Oxana Shevel of substantive responses that fall into the category of “Other” on V54. Those 11 respondents who responded “other” (V54=94) but did not name anything or anyone when asked what “other” meant (thus being missing values in V54othcod) are coded as “neutral.”*
generate odwhodiditx3=1 if V54==1
replace odwhodiditx3=1 if V54==2
replace odwhodiditx3=2 if V54==3
replace odwhodiditx3=2 if V54==4
replace odwhodiditx3=2 if V54==5
replace odwhodiditx3=3 if V54==95
replace odwhodiditx3=3 if V54==96
replace odwhodiditx3=3 if V54==97
replace odwhodiditx3=3 if V54==98
replace odwhodiditx3=V54othcod if V54==94
replace odwhodiditx3=3 if odwhodiditx3==.

*DV for logistic regression for aditional media analysis
generate odamdi=1 if odwhodiditx3==1
replace odamdi=0 if odamdi==.
 
generate odpmdi=1 if odwhodiditx3==2
replace odpmdi=0 if odpmdi==.

*Demographics
generate age05 = Vage_years
generate agegrp05 = Vage_6gr 
recode V58 (1=0) (2=1), gen(female)
generate commsize=Va14
recode V61 (1=1) (2=1) (3=2) (4=3) (5=4) (6=5) (7=6) (8=6) (97=.) (98=.), gen(educ05)

*Geographic place: Donbas, Galica, Odesa
recode Va13 (7=1) (13=1), gen(donbas)
replace donbas=0 if donbas~=1
recode Va13 (20=1) (11=1) (14=1), gen(galicia)
replace galicia=0 if galicia~=1
generate odesa=1 if Va13==16
replace odesa=0 if odesa==.

*Language of comfortable use
recode Vlang (1=0) (2=1) (3=0) (4=0) (5=0) (6=0), gen(rulangsvy05)
recode Vlang (1=1) (2=0) (3=0) (4=0) (5=0) (6=0), gen(uklangsvy05)

*Nationality
recode V64 (1=1) (2=0) (3=0) (7=0) (8=0) (.=0), gen(runats05)
recode V64 (1=0) (2=1) (3=0) (7=0) (8=0) (.=0), gen(uknats05)

*Orthodox Church - Moscow Patriarchate
generate orthmos05=1 if V73==2
replace orthmos05=0 if orthmos05==.

*Orthodox Church - Kyiv Patriarchate
generate orthkyiv05=1 if V73==1
replace orthkyiv05=0 if orthkyiv05==.

*Protest participation
generate maidan05=1 if V35_1==2 | V35_1==3 | V35_3==2 | V35_3==3
replace maidan05=0 if maidan05==.

generate antimaidan05=1 if V35_2==2 | V35_2==3 | V35_4==2 | V35_4==3 | V35_5==2 | V35_5==3
replace antimaidan05=0 if antimaidan05==.

*PARTISANSHIP: naming the party on question V13text (party member) or V16text or V19text (the usual my party series of questions). Respondents naming two parties are counted for neither party. Variables are created for the biggest parties, with at least 1% of respondents as partisans.
generate tpbat05=1 if V13text=="Batkivshchyna" | V16text=="Batkivshchyna" | V19text == "Batkivshchyna"
replace tpbat05=0 if tpbat05==. 
generate tppr05=1 if V13text=="Party of Regions" | V16text=="Party of Regions" | V19text == "Party of Regions"
replace tppr05=0 if tppr05==.

*Approval of Yatseniuk Performance as PM. This codes as 0 those who volunteer that he is not really/legitimately PM*
recode V20_1 (1=5) (2=4) (3=3) (4=2) (5=1) (6=0) (7=.) (8=.), gen(appyatspm05)

*Attitudes on regional language autonomy, EU, ATO
recode V57_5 (1=4) (2=3) (3=2) (4=1) (7=.) (8=.), gen(reglangaut05)
recode V48_4 (1=4) (2=3) (3=2) (4=1) (7=.) (8=.), gen(proeu05)
recode V57_7 (1=4) (2=3) (3=2) (4=1) (7=.) (8=.), gen(ato05)

*MEDIA*
generate uktv05=1 if V4_1==1 | V4_1==2 | V4_1==3 | V4_1==4 | V4_1==5 | V4_1==6
replace uktv05=0 if uktv05==.
generate rftv05=1 if V4_7==1
replace rftv05=0 if rftv05==.
generate notv05=1 if V1==2
replace notv05=0 if notv05==.
generate itv05=1 if V4_10==1
replace itv05=0 if itv05==.
generate odnoklas05=1 if V8_1==1
replace odnoklas05=0 if odnoklas05==.
generate vk05=1 if V8_2==1
replace vk05=0 if vk05==.
generate facebook05=1 if V8_3==1
replace facebook05=0 if facebook05==.

*Calculations for the frequency of watching Russian/Ukrainian TV used in robustness checks were made as follows. Note: this measure is not direct since only the frequency of TV new watching generally is recorded; the variables reflect the frequency of TV news watching among those who, respectively, say they watch Ukrainian and RussianTV.*
recode V3 (1=4) (2=3) (3=2) (4=1) (5=0), gen(tvfreq05)
replace tvfreq05=0 if V1==2
replace tvfreq05=. if V3==7
generate uktvfreq05=uktv05*tvfreq05
generate rftvfreq05=rftv05*tvfreq05


save "/Users/olaonuch/Desktop/UCEPS_ODESA_WITHVsRECODED_20171125.dta"

*Correlations matrix*

correlate rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 uklangsvy05 runats05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize


*Main regression

eststo clear

*DV: 1 = pro-Russians did it, 2 = pro-Ukrainians did it, 3 = other answers (baseline category)*

*Substituting means for missings for IVs*

qui: svy: mean (educ05)
qui: replace educ05=_b[educ05] if educ05==.
qui: svy: mean (appyatspm05)
qui: replace appyatspm05=_b[appyatspm05] if appyatspm05==.
qui: svy: mean (ato05)
qui: replace ato05=_b[ato05] if ato05==.
qui: svy: mean (proeu05)
qui: replace proeu05=_b[proeu05] if proeu05==.
qui: svy: mean (reglangaut05)
qui: replace reglangaut05=_b[reglangaut05] if reglangaut05==.

*Standardizing all IVs to 0-1 that are not already. The suffix ìuî on the var name stands for ìunit.î*

qui: summarize age05, meanonly
qui: gen age05u=(age05-r(min))/(r(max)-r(min))
qui: summarize educ05, meanonly
qui: gen educ05u=(educ05-r(min))/(r(max)-r(min))
qui: summarize commsize, meanonly
qui: gen commsizeu=(commsize-r(min))/(r(max)-r(min))
qui: summarize appyatspm05, meanonly
qui: gen appyatspm05u=(appyatspm05-r(min))/(r(max)-r(min))
qui: summarize ato05, meanonly
qui: gen ato05u=(ato05-r(min))/(r(max)-r(min))
qui: summarize proeu05, meanonly
qui: gen proeu05u=(proeu05-r(min))/(r(max)-r(min))
qui: summarize reglangaut05, meanonly
qui: gen reglangaut05u=(reglangaut05-r(min))/(r(max)-r(min))


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05u ato05u proeu05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(1)) post vce(unconditional)


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05u ato05u proeu05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(2)) post vce(unconditional)


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05u ato05u proeu05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(3)) post vce(unconditional)

esttab using ODtab320171021.rtf, replace compress b(%9.2f) star(* 0.05 ** 0.01) se wide varwidth(22) coeflabel(uktv05 "Ukrainian TV news" itv05 "Internet TV" rftv05 "Russian TV news" facebook05 "Facebook" vk05 "VKontakte" odnoklas05 "Odnoklassniki" notv05 "No TV" odesa "Odesa" tpbat05 "Batkivshchyna Party" tppr05 "Party of Regions" maidan05 "Maidan participant" antimaidan05 "Antimaidan participant" appyatspm05u "Approves Yatseniuk" ato05u "Pro-ATO" proeu05u "Pro-EU" reglangaut05u "Language autonomy" rulangsvy05 "Russian-speaker" runats05 "Russian ethnicity" orthmos05 "Orthodox (Moscow)" orthkyiv05 "Orthodox (Kyiv)" donbas "Donbas" galicia "Galicia" age05u "Age" educ05u "Education" female "Female" commsizeu "Larger community") nonumbers mtitles("Anti-Maidan Did It" "Pro-Maidan Did It" "Neutral") title("Table 3. Full effect of factors on probability of adopting beliefs about Odesa killings") addnote("Note: Multinomial logit.")

eststo clear


*Short regression = same as main regression except it drops values/attitudes & media

eststo clear

svy: mlogit odwhodiditx3 tpbat05 tppr05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(1)) post vce(unconditional)


svy: mlogit odwhodiditx3 tpbat05 tppr05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(2)) post vce(unconditional)


svy: mlogit odwhodiditx3 tpbat05 tppr05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(3)) post vce(unconditional)

esttab using ODtab3b20171021.rtf, replace compress b(%9.2f) star(* 0.05 ** 0.01) se wide varwidth(22) coeflabel(tpbat05 "Batkivshchyna Party" tppr05 "Party of Regions" rulangsvy05 "Russian-speaker" runats05 "Russian ethnicity" orthmos05 "Orthodox (Moscow)" orthkyiv05 "Orthodox (Kyiv)" donbas "Donbas" galicia "Galicia" age05u "Age" educ05u "Education" female "Female" commsizeu "Larger community") nonumbers mtitles("Anti-Maidan Did It" "Pro-Maidan Did It" "Neutral") title("Table 3b. Full effect of factors on probability of adopting beliefs about the Odesa killings") addnote("Note: Multinomial logit.")

eststo clear



*Creating table summarizing frequencies of variables of interest.*

use "/Users/olaonuch/Desktop/UCEPS_ODESA_WITHVsRECODED_20171125.dta", clear

la var rftv05 "Watches Russian TV news"
la def ny 0 "No" 1 "Yes", modify 
la val rftv05 ny

la var uktv05 "Watches Ukrainian TV news"
la val uktv05 ny

la var notv05 "No TV"
la def notv 0 "Watches TV" 1 "No TV", modify 
la val notv05 notv

la var itv05 "Watches Internet TV news"
la val itv05 ny

la var facebook05 "Facebook user"
la val facebook05 ny

la var vk05 "VKontakte user"
la val vk05 ny

la var odnoklas05 "Odnoklassniki user"
la val odnoklas05 ny

la var tpbat05 "Batkivshchyna partisan"
la val tpbat05 ny

la var tppr05 "Party of Regions partisan"
la val tppr05 ny

la var maidan05 "Maidan participant"
la val maidan05 ny

la var antimaidan05 "Antimaidan participant"
la val antimaidan05 ny

la var appyatspm05 "Approves Yatseniuk work as PM"
la def appyatspm 0 "Illegitimate" 1 "Absolutely not" 2 "Mostly not" 3 "Mixed" 4 "Mostly yes" 5 "Fully yes", modify 
la val appyatspm appyatspm

la var ato05 "Attitude to ATO"
la def agree4pt 1 "Fully disagree" 2 "Tend to disagree" 3 "Tend to agree" 4 "Fully agree", modify 
la val ato05 agree4pt

la var proeu05 "Attitude to EU"
la val proeu05 agree4pt

la var reglangaut05 "Attitude to regional language autonomy"
la val reglangaut05 agree4pt

la var Vlang "Language of Comfortable Use"
la def Vlang2 1 "Russian" 2 "Ukrainian" 3 "More Ukrainian" 4 "Both (answers Ukr.)" 5 "More Russian" 6 "Both (answers Rus.)", modify 
la val Vlang Vlang2

la var runats05 "Russian (mainly)"
la val runats05 ny

la var uknats05 "Ukrainian (mainly)"
la val uknats05 ny

la var orthmos05 "Orthodox (Moscow)"
la val orthmos05 ny

la var orthkyiv05 "Orthodox (Kyiv)"
la val orthkyiv05 ny

la var odesa "Odesa"
la val odesa ny

la var donbas "Donbas"
la val donbas ny

la var galicia "Galicia"
la val galicia ny

la var agegrp05 "Age group"
la def agegrp 1 "Under 30" 2 "30-39" 3 "40-49" 4 "50-59" 5 "60-69" 6 "Over 70"
la val agegrp05 agegrp

la var educ05 "Education level"
la def educ 1 "Elementary or less" 2 "Incomplete secondary" 3 "Secondary" 4 "Specialized secondary" 5 "Incomplete higher" 6 "Higher"
la val educ05 educ

la var female "Female"
la val female ny

la var commsize "Lives in larger community"
la def commsize 1 "Rural" 2 "Settlement" 3 "Towns under 20,000" 4 "Small city (20,000-49,999)" 5 "City (50,000-99,999)" 6 "Big city (100-499,999)" 7 "Very big city (500,000 or more)"
la val commsize commsize

tabout rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 Vlang runats05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia agegrp05 educ05 female commsize using ODfreqtab20171021.rtf, replace cells(col) oneway mi ptotal(none) svy percent


qui: svy: mean (educ05)
qui: replace educ05=_b[educ05] if educ05==.
qui: svy: mean (appyatspm05)
qui: replace appyatspm05=_b[appyatspm05] if appyatspm05==.
qui: svy: mean (ato05)
qui: replace ato05=_b[ato05] if ato05==.
qui: svy: mean (proeu05)
qui: replace proeu05=_b[proeu05] if proeu05==.
qui: svy: mean (reglangaut05)
qui: replace reglangaut05=_b[reglangaut05] if reglangaut05==.


*First, Ukrainian/Russian TV frequency and prevalence of AMDI/PMDI*

svy: logit odamdi rftv05 i.uktvfreq05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins uktvfreq05, vce(unconditional)

marginsplot, yline(0) ylabel(.8 "80%" .7 "70%" .6 "60%" .5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) xlabel(0 "None" 1 "Once" 2 "Several times" 3 "Almost daily" 4 "Daily") title("a. Ukrainian TV news viewers", size(large)) ytitle("% believing AMDI", size(large)) xtitle("Frequency in last 7 days", size(large)) recast(scatter) xscale(range(-.25 4.25)) graphr(color(white)) saving(ODuktvfreqamdi20171021, replace)

svy: logit odpmdi rftv05 i.uktvfreq05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins uktvfreq05, vce(unconditional)

marginsplot, yline(0) ylabel(.8 "80%" .7 "70%" .6 "60%" .5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) xlabel(0 "None" 1 "Once" 2 "Several times" 3 "Almost daily" 4 "Daily") title("b. Ukrainian TV news viewers", size(large)) ytitle("% believing PMDI", size(large)) xtitle("Frequency in last 7 days", size(large)) recast(scatter) xscale(range(-.25 4.25)) graphr(color(white)) saving(ODuktvfreqpmdi20171021, replace)


svy: logit odamdi i.rftvfreq05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins rftvfreq05, vce(unconditional)

marginsplot, yline(0) ylabel(.8 "80%" .7 "70%" .6 "60%" .5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) xlabel(0 "None" 1 "Once" 2 "Several times" 3 "Almost daily" 4 "Daily") title("c. Russian TV news viewers", size(large)) ytitle("% believing AMDI", size(large)) xtitle("Frequency in last 7 days", size(large)) recast(scatter) xscale(range(-.25 4.25)) graphr(color(white)) saving(ODrftvfreqamdi20171021, replace)


svy: logit odpmdi i.rftvfreq05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins rftvfreq05, vce(unconditional)

marginsplot, yline(0) ylabel(.8 "80%" .7 "70%" .6 "60%" .5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) xlabel(0 "None" 1 "Once" 2 "Several times" 3 "Almost daily" 4 "Daily") title("d. Russian TV news viewers", size(large)) ytitle("% believing PMDI", size(large)) xtitle("Frequency in last 7 days", size(large)) recast(scatter) xscale(range(-.25 4.25)) graphr(color(white)) saving(ODrftvfreqpmdi20171021, replace)


gr combine  ODuktvfreqamdi20171021.gph ODuktvfreqpmdi20171021.gph ODrftvfreqamdi20171021.gph ODrftvfreqpmdi20171021.gph, col(2) iscale(.5) title("Figure X. Estimated % believing AMDI & PMDI by frequency of watching TV news", size(medsmall)) graphr(color(white)) saving(ODfigtvfreq20171021, replace)



*Second, Effect of Ukrainian/Russian TV channels on AMDI/PMDI interacted with education.*

svy: logit odamdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("b. Ukrainian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukpmdi20171021, replace)


svy: logit odamdi rftv05##c.educ05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize

margins, dydx(rftv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvruamdi20171021, replace)


svy: logit odpmdi rftv05##c.educ05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize

margins, dydx(rftv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("d. Russian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvrupmdi20171021, replace)


gr combine  ODedxtvukamdi20171021.gph ODedxtvukpmdi20171021.gph ODedxtvruamdi20171021.gph ODedxtvrupmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by education", size(medsmall)) graphr(color(white)) saving(ODfigedxtv20171021, replace)





*SECOND-PLUS: Effect of Ukrainian TV channels on AMDI/PMDI interacted with education holding ethnicity constant.*

svy: logit odamdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional) at(runats05=1)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Effect on AMDI among Russians", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukxruamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional) at(runats05=1)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("b. Effect on PMDI among Russians", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukxrupmdi20171021, replace)


svy: logit odamdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional) at(runats05=0)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("c. Effect on AMDI among non-Russians", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukxukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.educ05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 female commsize
 
margins, dydx(uktv05) at(educ05=1 educ05=2 educ05=3 educ05=4 educ05=5 educ05=6) vce(unconditional) at(runats05=0)

marginsplot, yline(0) ylabel(.5 "50%" .4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("d. Effect on PMDI among non-Russians", size(large)) ytitle("Full effect", size(large)) xtitle("Education level", size(large)) xlabel(1 "Lowest" 6 "Highest") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODedxtvukxukpmdi20171021, replace)


gr combine  ODedxtvukxruamdi20171021.gph ODedxtvukxrupmdi20171021.gph ODedxtvukxukamdi20171021.gph ODedxtvukxukpmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of Ukrainian TV news on believing AMDI/PMDI" "by education and ethnicity", size(medsmall)) graphr(color(white)) saving(ODfigedxtvxid20171021, replace)



*Third, Ukrainian/Russian TV channels interacted with Russian/Ukrainian identity on AMDI/PMDI.*

svy: logit odamdi rftv05##runats05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(runats05=0 runats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Non-Russian" 1 "Russian") title("a. Russian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidruruamdi20171021, replace)


svy: logit odpmdi rftv05##runats05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(runats05=0 runats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Non-Russian" 1 "Russian") title("b. Russian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidrurupmdi20171021, replace)


svy: logit odamdi rftv05##uknats05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uklangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(uknats05=0 uknats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Non-Ukrainian" 1 "Ukrainian") title("c. Russian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidruukamdi20171021, replace)


svy: logit odpmdi rftv05##uknats05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uklangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(uknats05=0 uknats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Non-Ukrainian" 1 "Ukrainian") title("d. Russian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidruukpmdi20171021, replace)


svy: logit odamdi rftv05 uktv05##runats05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(runats05=0 runats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Non-Russian" 1 "Russian") title("e. Ukrainian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidukruamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##runats05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(runats05=0 runats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Non-Russian" 1 "Russian") title("f. Ukrainian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidukrupmdi20171021, replace)


svy: logit odamdi rftv05 uktv05##uknats05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uklangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(uknats05=0 uknats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Non-Ukrainian" 1 "Ukrainian") title("g. Ukrainian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidukukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##uknats05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uklangsvy05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(uknats05=0 uknats05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Non-Ukrainian" 1 "Ukrainian") title("h. Ukrainian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Self-stated ethnicity", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxidukukpmdi20171021, replace)


gr combine ODtvxidruruamdi20171021.gph  ODtvxidrurupmdi20171021.gph ODtvxidruukamdi20171021.gph ODtvxidruukpmdi20171021.gph ODtvxidukruamdi20171021.gph  ODtvxidukrupmdi20171021.gph  ODtvxidukukamdi20171021.gph  ODtvxidukukpmdi20171021.gph  , col(4) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by ethnicity", size(medsmall)) graphr(color(white)) saving(ODfigtvxid20171021, replace)



*Fourth, by TV effects by language use*

svy: logit odamdi rftv05##rulangsvy05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(rulangsvy05=0 rulangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Not Russian" 1 "Russian") title("a. Russian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangruruamdi20171021, replace)


svy: logit odpmdi rftv05##rulangsvy05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(rulangsvy05=0 rulangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Not Russian" 1 "Russian") title("b. Russian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangrurupmdi20171021, replace)


svy: logit odamdi rftv05##uklangsvy05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(uklangsvy05=0 uklangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Not Ukrainian" 1 "Ukrainian") title("c. Russian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangruukamdi20171021, replace)


svy: logit odpmdi rftv05##uklangsvy05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(rftv05) at(uklangsvy05=0 uklangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%" -.3 "-30%" -.4 "-40%" -.5 "-50%",angle(horizontal)) xlabel(0 "Not Ukrainian" 1 "Ukrainian") title("d. Russian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangruukpmdi20171021, replace)


svy: logit odamdi rftv05 uktv05##rulangsvy05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(rulangsvy05=0 rulangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Not Russian" 1 "Russian") title("e. Ukrainian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangukruamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##rulangsvy05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(rulangsvy05=0 rulangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Not Russian" 1 "Russian") title("f. Ukrainian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangukrupmdi20171021, replace)


svy: logit odamdi rftv05 uktv05##uklangsvy05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(uklangsvy05=0 uklangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Not Ukrainian" 1 "Ukrainian") title("g. Ukrainian TV on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangukukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##uklangsvy05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 uknats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(uktv05) at(uklangsvy05=0 uklangsvy05=1) vce(unconditional)

marginsplot, yline(0) ylabel(.3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) xlabel(0 "Not Ukrainian" 1 "Ukrainian") title("h. Ukrainian TV on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Language", size(large)) recast(scatter) xscale(range(-.25 1.25)) graphr(color(white)) saving(ODtvxlangukukpmdi20171021, replace)


gr combine ODtvxlangruruamdi20171021.gph  ODtvxlangrurupmdi20171021.gph ODtvxlangruukamdi20171021.gph ODtvxlangruukpmdi20171021.gph ODtvxlangukruamdi20171021.gph  ODtvxlangukrupmdi20171021.gph  ODtvxlangukukamdi20171021.gph  ODtvxlangukukpmdi20171021.gph  , col(4) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by language", size(medsmall)) graphr(color(white)) saving(ODfigtvxlang20171021, replace)



*Fifth, effect of Ukrainian/Russian TV channels on AMDI/PMDI interacted with age group (1-6).*

svy: logit odamdi rftv05 uktv05##c.agegrp05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia educ05 female commsize
 
margins, dydx(uktv05) at(agegrp05=1 agegrp05=2 agegrp05=3 agegrp05=4 agegrp05=5 agegrp05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Age group", size(large)) xlabel(1 "<30" 6 "70+") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODagextvukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.agegrp05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia educ05 female commsize
 
margins, dydx(uktv05) at(agegrp05=1 agegrp05=2 agegrp05=3 agegrp05=4 agegrp05=5 agegrp05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Age group", size(large)) xlabel(1 "<30" 6 "70+") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODagextvukpmdi20171021, replace)


svy: logit odamdi rftv05##c.agegrp05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia educ05 female commsize

margins, dydx(rftv05) at(agegrp05=1 agegrp05=2 agegrp05=3 agegrp05=4 agegrp05=5 agegrp05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Age group", size(large)) xlabel(1 "<30" 6 "70+") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODagextvruamdi20171021, replace)


svy: logit odpmdi rftv05##c.agegrp05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia educ05 female commsize

margins, dydx(rftv05) at(agegrp05=1 agegrp05=2 agegrp05=3 agegrp05=4 agegrp05=5 agegrp05=6) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Age group", size(large)) xlabel(1 "<30" 6 "70+") recast(scatter) xscale(range(.75 6.25)) graphr(color(white)) saving(ODagextvrupmdi20171021, replace)


gr combine  ODagextvukamdi20171021.gph ODagextvukpmdi20171021.gph ODagextvruamdi20171021.gph ODagextvrupmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by age", size(medsmall)) graphr(color(white)) saving(ODfigagextv20171021, replace)



*Sixth, effect of Ukrainian/Russian TV channels on AMDI/PMDI interacted with Approves Yatseniuk (strongest divider among attitudes).*

svy: logit odamdi rftv05 uktv05##c.appyatspm05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(appyatspm05=0 appyatspm05=1 appyatspm05=2 appyatspm05=3 appyatspm05=4 appyatspm05=5) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Approve Yatseniuk?", size(large)) xlabel(0 "Illegitimate" 5 "Fully Approve") recast(scatter) xscale(range(-.25 5.25)) graphr(color(white)) saving(ODyatsxtvukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.appyatspm05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(appyatspm05=0 appyatspm05=1 appyatspm05=2 appyatspm05=3 appyatspm05=4 appyatspm05=5) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Approve Yatseniuk?", size(large)) xlabel(0 "Illegitimate" 5 "Fully Approve") recast(scatter) xscale(range(-.25 5.25)) graphr(color(white)) saving(ODyatsxtvukpmdi20171021, replace)


svy: logit odamdi rftv05##c.appyatspm05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(appyatspm05=0 appyatspm05=1 appyatspm05=2 appyatspm05=3 appyatspm05=4 appyatspm05=5) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Approve Yatseniuk?", size(large)) xlabel(0 "Illegitimate" 5 "Fully Approve") recast(scatter) xscale(range(-.25 5.25)) graphr(color(white)) saving(ODyatsxtvruamdi20171021, replace)


svy: logit odpmdi rftv05##c.appyatspm05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05 proeu05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(appyatspm05=0 appyatspm05=1 appyatspm05=2 appyatspm05=3 appyatspm05=4 appyatspm05=5) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Approve Yatseniuk?", size(large)) xlabel(0 "Illegitimate" 5 "Fully Approve") recast(scatter) xscale(range(-.25 5.25)) graphr(color(white)) saving(ODyatsxtvrupmdi20171021, replace)


gr combine  ODyatsxtvukamdi20171021.gph ODyatsxtvukpmdi20171021.gph ODyatsxtvruamdi20171021.gph ODyatsxtvrupmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by attitude to Yatseniuk", size(medsmall)) graphr(color(white)) saving(ODfigyatsxtv20171021, replace)


*Seventh, effect of Ukrainian/Russian TV channels on AMDI/PMDI interacted with support for pro-EU attitudes*

svy: logit odamdi rftv05 uktv05##c.proeu05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(proeu05=1 proeu05=2 proeu05=3 proeu05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Join EU?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODeuxtvukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.proeu05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(proeu05=1 proeu05=2 proeu05=3 proeu05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Join EU?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODeuxtvukpmdi20171021, replace)


svy: logit odamdi rftv05##c.proeu05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(proeu05=1 proeu05=2 proeu05=3 proeu05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Join EU?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODeuxtvruamdi20171021, replace)


svy: logit odpmdi rftv05##c.proeu05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 reglangaut05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(proeu05=1 proeu05=2 proeu05=3 proeu05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Join EU?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODeuxtvrupmdi20171021, replace)


gr combine  ODeuxtvukamdi20171021.gph ODeuxtvukpmdi20171021.gph ODeuxtvruamdi20171021.gph ODeuxtvrupmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by EU support", size(medsmall)) graphr(color(white)) saving(ODfigeuxtv20171021, replace)




*Eighth, effect of Ukrainian/Russian TV channels on AMDI/PMDI interacted with support for pro-language autonomy attitudes*

svy: logit odamdi rftv05 uktv05##c.reglangaut05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(reglangaut05=1 reglangaut05=2 reglangaut05=3 reglangaut05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Regional language autonomy?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODautxtvukamdi20171021, replace)


svy: logit odpmdi rftv05 uktv05##c.reglangaut05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize
 
margins, dydx(uktv05) at(reglangaut05=1 reglangaut05=2 reglangaut05=3 reglangaut05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.4 "40%" .3 "30%" .2 "20%" .1 "10%" 0 "0" -.1 "-10%",angle(horizontal)) title("a. Ukrainian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Regional language autonomy?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODautxtvukpmdi20171021, replace)


svy: logit odamdi rftv05##c.reglangaut05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(reglangaut05=1 reglangaut05=2 reglangaut05=3 reglangaut05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on AMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Regional language autonomy?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODautxtvruamdi20171021, replace)


svy: logit odpmdi rftv05##c.reglangaut05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 appyatspm05 ato05 proeu05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05 educ05 female commsize

margins, dydx(rftv05) at(reglangaut05=1 reglangaut05=2 reglangaut05=3 reglangaut05=4) vce(unconditional)

marginsplot, yline(0) ylabel(.2 "20%" .1 "10%" 0 "0" -.1 "-10%" -.2 "-20%",angle(horizontal)) title("c. Russian TV news on PMDI", size(large)) ytitle("Full effect", size(large)) xtitle("Regional language autonomy?", size(large)) xlabel(1 "Fully disagree" 4 "Fully agree") recast(scatter) xscale(range(.75 4.25)) graphr(color(white)) saving(ODautxtvrupmdi20171021, replace)


gr combine  ODautxtvukamdi20171021.gph ODautxtvukpmdi20171021.gph ODautxtvruamdi20171021.gph ODautxtvrupmdi20171021.gph, col(2) iscale(.5) title("Figure X. Full effect of TV news on believing AMDI/PMDI by language autonomy view", size(medsmall)) graphr(color(white)) saving(ODfigautxtv20171021, replace)




*ROBUSTNESS CHECK: Now running the main regression, but substituting ukrainian language and ethnicity for russian language and ethicity*



*ROBUSTNESS CHECK: Now running the main regression, but dropping Approve Yatseniuk and Pro-EU to see if Batkivshchyna partisanship is significant in this case*

use "/Users/olaonuch/Desktop/UCEPS_ODESA_WITHVsRECODED_20171125.dta", clear

*DV: 1 = pro-Russians did it, 2 = pro-Ukrainians did it, 3 = other answers (baseline category)*

*Substituting means for missings for IVs*

qui: svy: mean (educ05)
qui: replace educ05=_b[educ05] if educ05==.
qui: svy: mean (appyatspm05)
qui: replace appyatspm05=_b[appyatspm05] if appyatspm05==.
qui: svy: mean (ato05)
qui: replace ato05=_b[ato05] if ato05==.
qui: svy: mean (proeu05)
qui: replace proeu05=_b[proeu05] if proeu05==.
qui: svy: mean (reglangaut05)
qui: replace reglangaut05=_b[reglangaut05] if reglangaut05==.

*Standardizing all IVs to 0-1 that are not already. The suffix ìuî on the var name stands for ìunit.î*

qui: summarize age05, meanonly
qui: gen age05u=(age05-r(min))/(r(max)-r(min))
qui: summarize educ05, meanonly
qui: gen educ05u=(educ05-r(min))/(r(max)-r(min))
qui: summarize commsize, meanonly
qui: gen commsizeu=(commsize-r(min))/(r(max)-r(min))
qui: summarize appyatspm05, meanonly
qui: gen appyatspm05u=(appyatspm05-r(min))/(r(max)-r(min))
qui: summarize ato05, meanonly
qui: gen ato05u=(ato05-r(min))/(r(max)-r(min))
qui: summarize proeu05, meanonly
qui: gen proeu05u=(proeu05-r(min))/(r(max)-r(min))
qui: summarize reglangaut05, meanonly
qui: gen reglangaut05u=(reglangaut05-r(min))/(r(max)-r(min))


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(1)) post vce(unconditional)


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(2)) post vce(unconditional)


svy: mlogit odwhodiditx3 rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05 tpbat05 tppr05 maidan05 antimaidan05 ato05u reglangaut05u rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu, base(3)

eststo: margins, dydx(*) predict(outcome(3)) post vce(unconditional)


esttab using ODtab3dropeuyats20171021.rtf, replace compress b(%9.2f) star(* 0.05 ** 0.01) se wide varwidth(22) coeflabel(uktv05 "Ukrainian TV news" itv05 "Internet TV" rftv05 "Russian TV news" facebook05 "Facebook" vk05 "VKontakte" odnoklas05 "Odnoklassniki" notv05 "No TV" odesa "Odesa" tpbat05 "Batkivshchyna Party" tppr05 "Party of Regions" maidan05 "Maidan participant" antimaidan05 "Antimaidan participant" ato05u "Pro-ATO" reglangaut05u "Language autonomy" rulangsvy05 "Russian-speaker" runats05 "Russian ethnicity" orthmos05 "Orthodox (Moscow)" orthkyiv05 "Orthodox (Kyiv)" donbas "Donbas" galicia "Galicia" age05u "Age" educ05u "Education" female "Female" commsizeu "Larger community") nonumbers mtitles("Anti-Maidan Did It" "Pro-Maidan Did It" "Neutral") title("Table 3. Full effect of factors on probability of adopting beliefs about the May 2014 Odesa killings") addnote("Note: Multinomial logit.")

eststo clear



*MEDIA AS DVs (rftv05 uktv05 notv05 itv05 facebook05 vk05 odnoklas05) with demographic IVs.*

svy: logit rftv05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit uktv05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit notv05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit itv05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit facebook05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit vk05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

svy: logit odnoklas05 rulangsvy05 runats05 orthmos05 orthkyiv05 odesa donbas galicia age05u educ05u female commsizeu

eststo: margins, dydx(*) post vce(unconditional)

esttab using ODmediadv20171021.rtf, replace compress b(%9.2f) star(* 0.05 ** 0.01) se varwidth(22) coeflabel(odesa "Odesa" rulangsvy05 "Russian-speaker" runats05 "Russian ethnicity" orthmos05 "Orthodox (Moscow)" orthkyiv05 "Orthodox (Kyiv)" donbas "Donbas" galicia "Galicia" age05u "Age" educ05u "Education" female "Female" commsizeu "Larger community") nonumbers mtitles("RuTV" "UkTV" "NoTV" "ITV" "FB" "VK" "Odn") title("Table X. Full effect of demographics on probability of consuming different media May 2014") addnote("Note: logit.")

eststo clear




log close
